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We present a self consistent approach to Coulomb gauge Hamiltonian QCD 



> 

m 
oo 
o 
in 
o 
o 
o 



which allows one to relate single gluon spectral properties to the long range 



behavior of the confining interaction. Nonperturbative renormalization is dis- 



cussed. The numerical results are in good agreement with phenomenological 




and lattice forms of the static potential. 
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I. INTRODUCTION 



There is a long history |T]-[§] of using mean field techniques to explore dynamical chiral 
symmetry breaking in Quantum Chromodynamics. Chiral symmetry breaking is of funda- 
mental importance because it is directly related to the vacuum structure of QCD. It also 
has important practical consequences, such as making chiral perturbation theory a viable 
approximation to low energy QCD. Mean field equations are typically derived by analyzing 
a truncation of the Schwinger- Dyson equations or, when dealing with a Hamiltonian, by 
performing a Bogoliubov-Valatin transformation on the particle basis. 

An added complexity with Hamiltonian approaches has to do with renormalization. Un- 
like covariant approaches to field theory which generate a finite number of counterterms, 
the Hamiltonian formalism necessarily involves a noncovariant truncation of the theory and 
hence may generate noncanonical counterterms. As a result it may not be possible to remove 
the ultraviolet cutoff. In this sense, an effective field theory approach to Hamiltonian-based 
renormalization becomes natural. 

The problem of renormalization in models derived from a QCD Hamiltonian is often 
ignored. An exception was the work of Adler and Davis [1]] who have noted that, when 
using the BCS Ansatz, consistency with the Ward identities imposes a definite counterterm 
structure on the gap equation. They also noted that calculations should be performed with 
the full Hamiltonian (as opposed to the normal ordered Hamiltonian). 

Although the quark sector of the vacuum has been studied relatively well, the gluonic 
sector has largely been ignored. However, the gluonic vacuum is also relevant to studies of 
chiral symmetry breaking because the quark and gluonic gap equations are coupled and, more 
importantly, the gluon propagator is related to the quark- ant iquark confining interaction. 
Furthermore, the gluonic vacuum determines the properties of the gluonic quasiparticles 
which in turn may be used to construct a Fock space expansion of hadrons such as glueballs 
or hybrids ||. 

In this paper we derive and renormalize the gap equation for the purely gluonic sector 
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of QCD. This work extends three previous papers. The first of these presented the unrenor- 
malized gluonic gap equation and calculated the gluon condensate and glueball spectrum 
||. The other two |||7| described how to renormalize Hamiltonian QCD using the similarity 
or flow equation methods for evolving the renormalization scale. Both methods are neces- 
sarily perturbative, which means that gap equations may be derived with one (or more) loop 
corrections. However, since the phenomena being studied are nonperturbative, it is prefer- 
able to construct a nonperturbative renormalization scheme from the start. Unfortunately, 
this means abandoning the elegant flow equation methodology which is so well suited to 
Hamiltonian-based perturbation theory. 

We further extend the previous work by showing how to obtain the long range interaction 
between color sources. This is done in the BCS Ansatz by simultaneously considering the 
single gluon spectral function and the non-Abelian Coulomb potential. This subject was 
previously considered in Ref. and || (we discuss the various approaches below). We 
derive a self consistent gap equation for the gluon spectral function. Self consistency arises 
because the non-Abelian potential appearing in this equation depends on the gluon spectral 
function. The solution to the gap equation gives rise to an infrared enhancement in the 
effective non-Abelian potential which may be identified with linear confinement. 

In the following section we give a brief review of QCD in Coulomb gauge and present our 
regularization scheme. We then discuss the renormalization scheme adopted, the gap equa- 
tion, and the effective interaction emerging from the self consistent solution. We conclude 
and discuss future applications in Section III. 

II. QCD IN COULOMB GAUGE 

Understanding the properties of soft gluons is one of the major challenges of hadronic 
physics. It is natural to study soft gluons using the Hamiltonian formulation of QCD in a 
physical gauge such as the Coulomb gauge. This has many advantages: it is closest in spirit 
to quantum mechanical models of QCD, all degrees of freedom are physical, no additional 
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constraints are need on the Fock space, the norm is positive definite, spurious retardation 
effects are minimized, and confinement may be rigorously identified with the non-Abelian 
Coulomb interaction in the heavy quark limit jllj . Finally this is also the natural choice for 
the study of nonrelativistic bound states of a few constituent degrees of freedom. It is of 
relevance in the gluonic sector since one would expect the gluonic quasiparticles to be heavy 
due to strong confining interactions. 

The pure gauge QCD Hamiltonian in Coulomb gauge may be written as [11| H = 
H g + He, where the terms are given by, 

H g = Tr J d 3 x [j~ x U JU + B b] (1) 

He = \g 2 J d*x£yJ- l p a (?L)K a \-K,y)Jp\y) (2) 

The second term is the instantaneous non-Abelian Coulomb interaction for QCD. The kernel 
K is given by 

K ab ( X ,y) = (x,a|(V • D) _1 (— V 2 )(V • D)" 1 ^*) (3) 

and the color charge density p a is given by 

p a (x) = f abc A b {x) ■ IT(x). (4) 

When expanded in powers of the strong coupling, the Coulomb kernel contains infinitely 
many terms arising from the inverse of the adjoint covariant derivative: 

B ab = S ab v _ g fabc A c . y ^ 

At lowest order in the coupling the kernel is given by K a6 (x, y) = 5 ab / (47r|x — y|) as in QED. 
The Faddeev-Popov determinant appears in the kinetic energy and Coulomb interaction 
terms due to the curvature of the gauge manifold and is given by J = det [V • D]. Lastly, 
the components of the non-Abelian magnetic field are given by 

B? = e ijk (VjA% + 9 -f ahc A)Ai) . (6) 



It should be stressed that this is a bare Hamiltonian. Regularization leads to counterterms 
which in practice will depend on the approximation schemes used to diagonalize the Hamil- 
tonian. As stated above, our goal is to study the QCD vacuum with the aid of the BCS 
Ansatz. In the past we have chosen to regulate the Hamiltonian by restricting its matrix 
elements in the basis of eigenstates of the bare kinetic energy to a band-diagonal form ||[7| . 
Although an elegant formulation of the renormalization procedure is possible with this choice 
of basis, it gives rise to lengthy expressions because all calculations must be performed in 
the particle basis. Furthermore, such a scheme does not seem to be particularly relevant 
when dealing with nonperturbative renormalization since there is no particular advantage 
in using the bare basis. Hence, it is preferable to employ the field basis and we introduce 
a field-based regulator which is analogous to Schwinger's point splitting. This consists of 
smearing field operators over a small spatial region: 

Here A is the UV cutoff and the fields are effectively smeared over a distance 0(1/ A). The 
mode expansions are 

A ' (x) ^/(S7i!( a - (k)+a ' ( - k) ) e "* k ^ (8) 

ASM = / 0<^(<*(k) - «T(-k)) e-*"& (9) 

where ui (k) = |k| in the perturbative vacuum. We shall subsequently drop the tildes. 
Contractions of the field operators which are needed below are given as follows 

t bf v \\ _ tab f Aj(fc) Jk-(x-y)- fc2 



ww^iwiJr ' (10) 



d * k n r^ g°W o*-(*-y)-T 
(2tt 



<n?(x)Ii;(y)> = 5 ab I ^D y (*)=Hpe*^-S» (11) 



(A?(x)Il*(y)> = i5<* j ^ £<M e *><*-y)-& (12) 
where Dij(k) = — kikj is the transverse delta function associated with Coulomb gauge. 
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A. Nonperturbative Renormalization 



In order to reduce the cutoff dependence induced by regularization we allow the cou- 
plings to be A-dependent and add an infinite set of counterterms to the Hamiltonian. The 
counterterms are organized in powers of the cutoff, 

5H= ^T^ W - (13) 

n=-2 A™ 

As long at the cutoff only affects the operator products at short relative distances, the 
counterterms O must be local. Thus they may be classified according to their canonical 
dimension, n + 1. Furthermore they have to preserve the unbroken symmetries i.e. be 
rotationally invariant, color, spin, and flavor singlets. Notice that the series starts at or- 
der A 2 . This is because the lowest dimension operator satisfying the above restrictions is 
/ <ixA 2 (x) which is of dimension —1, (n— —2). The needed marginal and relevant operators 
contributing to the gluon sector of the Hamiltonian are 

SH= M^J d3xA2 + c J ^lj rf , xAv2A 

+ z - 1 w- l j* x iP + mzij <exB > (14) 

We have established contact with more traditional notation by defining M 2 (A) = A 2 c_ 2 
where c_ 2 is dimensionless. Operators of dimension three and five do not occur. Thus the 
next set of relevant operators are of dimension six: f abc f cde A ia U tb A jd W e , A ■ UA ■ II, AV A A, 
nV 2 II, f abc \l m A b ■ VIP C , combinations of six gluon fields and gradients, and so on. 

It is important to recall that in an effective field theory the cutoff is not removed. Rather 
the various coefficients are determined by requiring that observables are accurate to a given 
order in p/A, where p is some characteristic, measurable momentum scale. In our case the 
c n (A) should be tuned to reproduce experimental data calculated in some nonperturbative 
scheme. Here we are referring to a number of standard many-body calculational schemes 
such as the BCS vacuum Ansatz, the Tamm-Dancoff truncation, or the random phase ap- 
proximation. 
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Transverse gluons do not contribute to the counterterms because we employ a BCS 
Ansatz in this work. Thus at this stage of the calculation, the A dependence of the coun- 
terterms will not match that of perturbative QCD. Nevertheless couplings to the transverse 
gluons with momenta below the cutoff are still present in the Hamiltonian and lead to 
important effects which (if the valence sector dominates) could be taken into account in per- 
turbation theory around quasiparticle bound states. Removal of the cutoff without taking 
into account the effects from transverse gluons would clearly be incorrect. This will be seen 
explicitly below when we compare, in the weak coupling limit, the counterterms calculated 
with the BCS Ansatz with the results coming from perturbation theory. 



B. The Gap Equation 



There are many equivalent formulations of the BCS approach. The one appropriate to 
many-body physics is the Bogoliubov-Valatin (BV) canonical transformation on the particle 
operators. In our case this may be written as 



a\{k) = c{k)a\{k) + s(k)o%\-k) 
a?(k) = c(k)a?(k) + s(k)at(-k) 



(15) 
(16) 



where the rotation is parameterized in terms of an unknown gap function, u, as 

c(k) 



1 
2 



u (k) 



+ 



\ u{k) \ uj Q {k) 



s{k) 



1 

2 



u (k) 



(17) 



(18) 



^ u(k) \| uj (k) 

Because we have employed a field-based regulator, the effect of the BV transformation is 
to simply replace u)o(k) = k with ujq = oj(k) in the mode expansions of Eq. (^|). Thus all 
matrix elements evaluated in the BCS vacuum are simply related to those evaluated in the 
perturbative vacuum with the same replacement. 

The gap equation may now be obtained from the one-body portion of the QCD Hamil- 
tonian. Normal ordering with respect to the BCS vacuum yields 



H 



d 3 q 



k- 



E(q)a? (q)aj(q) + G(q) (a*(q)«*(-q) + ff.c. 



(19) 



where 



£(g) 



2w(g) 



+ Zg 2 + M 2 (A) + g 2 c (A) 



4 J (2vr) 3 



+ ^ / e~A* V(k + q) (1 + {k ■ qf \ 



u 2 {k)+uo\q) 
Lu(k) 



+ira s (A)N c 



d 3 k _V (3- (k-q) 2 ) 



(2tt)* 



u(k) 



(20) 



and 



G(q) 



Au{k) 



co 2 (q) 



+ Zg 2 + M 2 (A) + c (A)g 2 



Nc f d 3 k , 

+ T7 (2^)3 



" ' " " V(k + q)e-^ (l + (k-qf 



uo 2 {k) - uo 2 {q) 
uj(k) 



+ na s (A)N c 



d 3 k (3 - (k ■ q) 2 ) 



e a 



(2tt) 3 u{k) 
The spectral function u(k) is obtained from solving the gap equation 

G(q) = 



(21) 



(22) 



which demands that the quasiparticle, BCS vacuum decouples from states with pairs of 
gluons -a feature reminiscent of the constituent quark model. 

The static potential in these expressions is the expectation value of the non-Abelian 
Coulomb kernel in the BCS vacuum state, 



v(<D 



d 3 re- lCir (BCS\K ab (r, 0)\BCS). 



(23) 



At lowest order in the coupling and in the perturbative vacuum this potential is simply 
47ra s (A)/q 2 . It is possible to show that the non-Abelian Coulomb term describes the com- 



plete quark-antiquark interaction in the heavy quark limit ||10|| . Thus one may use lattice 
data to determine that the potential takes on the familiar Coulomb+linear form. It is nat- 
ural to assume that, with the appropriate color structure, this form holds for gluonic color 



sources as well (in fact this is confirmed by lattice calculations |15|). Thus, in previous 



calculations [K]]7]] we have simply replaced the expectation of the non-Abelian kernel with a 
linear potential. However, one of our current goals is to demonstrate that the BCS Ansatz 
and linear confinement can be described consistently by the same formalism. The central 
idea is that the gap equation contains a kernel which is itself the expectation value of the 
non-Abelian Coulomb interaction in the BCS vacuum. Thus the kernel is functionally de- 
pendent on the gap function and it is possible to obtain both the gap function and the 
effective potential. Since it is known that the non-Abelian Coulomb interaction gives rise to 
confinement in the heavy quark limit, one can hope that this procedure will yield a linear 
potential at large distances |§, @. 

In an attempt to gain some insight into this hypothesis, we consider perturbative cor- 
rections to the effective potential. As stated above, the non-Abelian Coulomb kernel is an 
operator which depends on the transverse degrees of freedom of the gauge field (cf. Eq. (|3])). 
It can be shown that the (3 function can be obtained from loop corrections to this ker- 
nel [[HJ. One finds that at 0(a s ), there are three contributions to (3 = —P /(2tt) ■ a 2 s \ 
Po — Pc + P g + P q = +4A" C + (—1) + (— 2nf/3). The first term, /3 C comes directly from 
the expansion of the non-Abelian Coulomb potential. The second and third terms are due 
to mixing with two-gluon and two-quark intermediate states respectively. In the following 
we consider the perturbative corrections to the Coulomb kernel since it gives rise to the 
majority of the pure gauge P function. 

The contribution from the transverse gluons in the expansion of the Coulomb kernel may 
be computed with the expansion 



V D v V D J 



S ab + 2 9l ^f abc A c • V + 3g 2 ^f acd A d • V^/ cfee A e • V + 

(24) 



with the cutoff dependence of the couplings, gi = g(A)Zj(A) to be determined. To order gf, 
the Coulomb interaction is given by 



Ana s 
Vc{q) = — 2- 



+ 16^ J (^ A ^o) 



(25) 



where a s (A) = g 2 (A)/Air and 

I(q, A; u ) = f ^ {1 ~^f )e ~T a - ( 26 ) 

In the BCS vacuum the same expression arises except that ojq is replaced with the un- 
known function, uj. Thus, as stated above, the solution to the gap equation depends on a 
potential which itself is a function of the gluon dispersion relation. It therefore becomes pos- 
sible to obtain the gluon gap function and the nonperturbative potential in a self consistent 
fashion. 

We pursue this scenario by defining a nonperturbative coupling via the relation 

VM) - ^f^- (27) 

Furthermore, we take part of the rainbow-ladder approximation to V c by summing all dia- 
grams proportional to I n . This leads to an expression familiar from the perturbative leading 
log approximation, 

Although this potential may be inserted directly into the gap equation we have found it 
useful to approximate it with the aid of the following substitution 



47T 

IM -»• -3- log 



A 2 



Mi) 2 . 

As shown in Fig. 1, this form is accurate to roughly 10% for q/A < 0.1. 



(29) 
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FIG. 1. 

(/M - log 



Relative error due to approximation defined in Eq. 



The ratio 



)/I[uj] is plotted as a function of q/A for to(q) 2 = M 2 + q 2 . The two curves 



correspond to M/A = 0.01 and M/A = 0.1 (larger difference for smaller q). 

With this expression for V c , the gap equation is given by 

u(q) 2 „ o , . 9/AN , k . 9 , M „ / ft fc 2 (3 - (fc • g) 2 ) 

-f- = V + Af (A) + c„(A)^ + ™,(A)JV C J ^ ±jJU 

N c r d 3 k _2£ T ~ r/1 « / ? , 2 \ uo 2 (k) - tu 2 (q) 



(30) 



Many of the terms of higher order in 1/A also contribute to the gap equation. Additional 
terms include c 2 {A)q 2 uj{q) 2 /A 2 , d 2 {A)q 4 /A 2 , and e 2 (A)/A 2 • N c /3 • / exp(-A; 2 /A 2 )(cu(A;) 2 - 
u(q) 2 )/u(k). As expected, these only affect the solution at short range, thus the long range 
behavior of the effective potential is completely specified by the nonperturbative model (the 
BCS Ansatz in this case). Nevertheless, one may get an indication about the viability of the 
low energy model by studying the coefficients, C2, e^, etc - a strong cutoff dependence at low 
momenta would indicate that the nonperturbative model is inadequate and that additional 
counterterms are therefore required. 

At this point the standard procedure would be to choose values for the cutoff and coeffi- 
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cients, solve the gap equation to obtain u>, use the solution to obtain, say, glueball massesQ, 
and vary the coefficients so that the predictions agree with experimental (or lattice) data. 
However, in this study we choose to perform a simpler analysis and fit the derived form of 
the effective potential, V e /f = V(u>) to the Wilson loop lattice potential. This is similar in 
spirit to a lattice renormalization procedure advocated by Lepage and Mackenzie [IB|. Di- 



mension six operators are neglected so that the couplings to be varied are the five coefficients 
Z, c_2, Co, Z 2 , and a s . First we note that not all of the coefficients are independent since 
one requires that the effective potential, and therefore u, is A independent (to leading order 
in p/A). Multiplying the gap equation by Z and introducing d s = Za s , c_2 = Z(c_2 + L), 
c = Z 2 + Zcq leads to 

u\q) = ~c (A)q 2 + A 2 MA) + j^y 2 f dkk 2 d(q ■ k)V c (k + q) (l + (q ■ k) 2 ) 

(31) 



where 



V M = — ^TaT^ — T^Y ( 32 ) 



and L is defined to be the constant, A-dependent contribution from the third term on the 
right hand side of Eq. ( |3T)D (this term comes from the A A operator in the Hamiltonian) . In 
obtaining Eq. ([H|) we have further chosen Z = Z\ which is required to make the effective 
interaction cutoff independent. 

We now define a QCD scale, /i, as 



4tt 



fi 2 = A 2 e (33) 

and set m 2 (q) = oo 2 (q) — q 2 . The expression for the running coupling becomes cutoff inde- 
pendent (as long as m(g) is A- independent) and is given by 



1 This procedure is followed, without renormalization, in Ref. 
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« e// (g 2 ) = ^r^, VM = 4TOe/ / (g2) ■ (34) 

We note that if m(0) = fi the long range behavior of the effective potential corresponds to 
a linear potential in position space. Equating the coefficient of the leading 1/g 4 behavior 
with the standard form 6nb/q 4 (b is the string tension), yields /i 2 = 3(1 + m')(3 c b/(8iT) where 
ml = dm 2 (0) / dq 2 . The two remaining coefficients in the gap equations can thus be fixed by 
demanding that the effective potential represents linear confinement with the right slope. In 
numerical computations, instead of fixing m' we have chosen to fix m 2 (q) at a finite value of 
g 2 , typically q 2 = 1 GeV 2 , which is numerically easier. The two renormalization conditions, 
one for m(0) and one for m! (or in general for m 2 (qo) and m 2 (gi)), can be implemented 
directly into the gap equation by two subtractions, 



g 2 -g 2 A> d/eV ^ w (jfc) + g 2 -g 2 yo d/eV ^'^ 



where cjq = ^(?o), = ^((Zi) and 



AT r 1 - fc^ ~ 
V(*, q) = t^- 2 J d(q ■ k)k 2 e-7? V c (k + q)(l + (q ■ k) 2 ) (36) 



(An) 2 J-i 

The original couplings may be reconstructed as follows: 



2 q&Z-qfcl q\ [°° u 2 (k)-ul ql [™ u 2 (k) - u 2 



A^_ 2 = ^~ u 2 T - dkV(k, g ) - / » \ + dkV(k, q x 



q\ -ql qf-qVo " g? - g^ ' 

(37) 



The single gluon energy is given by 

V(fc,g) 



£(g) = u(q) 



7-1 



dk- 



(39) 



u(k) 

If a linear potential is to be recovered (as it must because we are fitting to lattice data), 
one must have u = \i at g = 0, (see Eq. (0)). Thus the QCD scale \i must be related to the 
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string tension. As discussed above, instead of taking the limit qi — > q and fixing du 2 (0)/dq 2 
to reproduce the magnitude of string tension we keep qi finite, choose q± = 1 GeV and fit lo\ 
to reproduce the strength of the confining potential at low momentum. In our computations 
the one dimensional nonlinear integral (after numerical evaluation of the angular integrals) 
equation was solved numerically using a modified Levenberg-Marquardt algorithm. We 
note that the equations are extremely sensitive to the proper treatment of the IR singularity 
(which, however, is integrable) and that some solutions which have appeared in the literature 
are misleading. The results for a e ff(q) = q 2 V e f f(q) j 'An are shown in Fig. 2. We find that 
m' ~ (typically m! = 0.01) thus the QCD scale /i can be directly calculated from the 
physical string tension. Taking b = 0.18 GeV 2 yields fi 2 = 0.2579 GeV 2 . Finally, for large 
momenta the running of the effective coupling is given by 



47T 



<//(?) = —7 q2 y (40) 

Pc lOg ^ 2579Ge y2 ) 

Over the a range of momenta the full solution to the effective coupling (or potential) can 
be well approximated by a running coupling version of the Cornell potential 

/ \ 4tt 3 b a%(A) 3 b 

a efM) = — — —w + n~ = t ^mi^_ + o- ( 41 ) 









Pc log 


"<zi 4 





2*» ( 1 + » logE ^) 2,f 



Here jiQ is a parameter which allows interpolation between the infrared and ultraviolet 
limits of the potential. We find (fio/fi) 2 ~ 10. Essentially identical results are obtained 
upon comparison with the standard Coulomb+linear form of the Cornell potential. 
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FIG. 2. The effective potential derived from the self-consistent gap equation. The dashed 
line is the numerical solution to the gap equation with A = 10 GeV (A dependence is negligible 
and would not be visible on this scale). The solid line is the fit to the running Cornell potential 
as discussed in the text. The discrepancy at low momenta is due to the finite mesh size in the 
numerical calculations. 

That the form of the derived effective potential so nearly reproduces linear confinement 
at long range is a tantalizing indication that our goal of deriving the confinement potential 
in a self-consistent manner from the BCS vacuum Ansatz may be achieved. 

The derived gluon dispersion relation u>(q) is shown in Fig. 3. The doubly subtracted 
form of the gap equation removes A 2 and log(A) dependence from the gap equation and 
leaves sub leading terms of order q 2 / A 2 . This is reflected in the cutoff independence evident 
in Fig. 3 at low momentum. In Fig. 3 we also compare the full solution of the gap equation 
for A = 10 with the approximate one obtained using the "Coulomb+linear" potential with 
the parameters fitted to the self consistent solution as discussed above. Both of them lead to 
m' ~ and both are renormalized at the same point, qi = 1 GeV, however the intermediate 
q dependence is somewhat different. As seen from Fig. 2, this small difference does not alter 
the behavior of the effective interaction. 
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FIG. 3. m 2 (q) = uj 2 (q) — q 2 is a function of momentum. The dashed and dotted lines 
correspond to A = 10 GeV and A = 100 GeV respectively. The solid line is the Cornell potential 
fit to the full solution with A = 10 GeV. The dash-dotted line is m 2 (q) = q 2 . 



One may expect that u(q) should approach the free solution, u (q) = q at large momenta. 
However, this need not be true, both in perturbation theory, where logarithms can ruin 
the asymptotic behavior, and nonperturbatively. This is because the counterterms Z B 2 
and coAVA give rise to the term c$q 2 in the gap equation. Renormalizing cu(q) at a low 
momentum scales (fa, qi « A will result in a logarithmic dependence on A in c - Thus the 
effective gluon mass m 2 (q) = uj 2 (q) — q 2 will remain proportional to q 2 at large momenta. 

As mentioned above, this behavior is also true in perturbation theory. In the following 
we examine the small coupling limit of the counterterm coefficients to establish contact with 
perturbative QCD. Evaluating Eqns 37 and 38 to O (a s = ci s (A)) yields 



C-2 = 



ajNc 
An 



4 8 
3 ~ 3 



a s /3 



N, 



c 



(42) 



where the first contribution comes from the vacuum expectation value of the Coulomb 
potential and the second from the A 4 term in the Hamiltonian, and 
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Co = 1 + ^ A4±4 In * U ^ A In *L U ^ A In A 2 + finite. (43) 
4vr 15g 2 -g 2 UJ 4vr 15 U ?J 4vr 15 



In perturbation theory, to 0(a s ), A-dependence of Z and Co follows from Eq. (^) and the 
requirement that the gluon energy (an observable in perturbation theory) is independent of 
A. Thus from Eq. ( |3SD it follows that 



E(q) = q 



„ i a a JV c 4, A 2 
4vr 3 q 2 



(44) 



which leads to 



and finally 



^-lnA 2 + finite (45) 

4:71 3 



a s N c 32 2 

C ° = 4~7{~15 (46) 



These expressions agree with the corresponding expressions in a full perturbative Hamil- 
tonian QCD calculation [0. Note, however, that in Ref. [0 it has been shown that Cq = to 
0(a s ). This arises because the self energy contribution from a gluon loop with transverse 
gluon exchange cancels the contribution from the Coulomb and the A 4 vacuum expectation 
valuesQ. The BCS calculation only contains the latter two terms. In a nonperturbative 
calculation based on the BCS Ansatz transverse gluons are confined and the single gluon 
energy of Eq. ( j39D becomes IR singular. As a result, transverse gluons only affect observables 
in bound state calculations of color singlet states. 

In Fig. 4 we show the behavior of the counterterms, c_2 and c"o for small couplings. The 
points represent the full expressions obtained from Eqns 37 and 38. The lines are the 
perturbative results given above with the finite portions fit to the full results. As expected 
the agreement is very good for small coupling. Fig. 4a shows a dramatic dependence of 
— c~2 on the value of the UV cutoff. This dependence originates with the subleading 1/q 4 



Note that the UV cutoff A used in Ref. QJ is by a factor of v2 larger then the one used here. 
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behavior of the effective potential which makes a log(A) contribution to c_2- The solid lines 
in Fig. 4a incorporate this behavior. The logarithmic term dominates in the small coupling 
limit and leads to the large relative splitting seen in the figure. 




Na/47i NctAta 



FIG. 4. Counterterms as a function a s in the weak coupling limit. The two set of points 
correspond to A = 10 GeV (squares) and A = 100 GeV (triangles). Solid lines are the results of 
the leading order 0(a s ) perturbative approximation given by Eqs. @) and (|i^). 



C. Relationship to Other Approaches 



The Schwinger-Dyson formalism if often used to discuss the issues presented here [jT4| . pT9 
Advantages are that it is covariant and that calculations beyond the rainbow approximation 
are possible. However, the resulting integral equations are four dimensional and are therefore 
difficult to solve. They are also normally evaluated in Euclidean space which introduces 
interpretational difficulties, especially when dealing with confinement. Furthermore, the 
connection to the parton model is unclear since it is difficult to continue wavefunctions 
to Minkowski space. In the Hamiltonian framework color nonsinglets decouple from the 
physical spectrum as in the Schwinger-Dyson approach, but hadronic wavefunctions are 
well-defined, and thus a parton picture of high energy QCD is possible. 

We also note that mean field formalisms of Coulomb gauge QCD have been extensively 
studied in Refs. || and M. The numerical solution presented in Ref. M are however doubtful, 
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as they do not have the proper large momentum behavior. In Ref. the Coulomb potential 
was studied with the aid of a model of the gluon dispersion function (the gap equation was 
not solved). 

A calculation with similar goals to those expressed here has recently appeared fll~T 



This paper employs the flow equation renormalization methodology to derive a gluonic gap 
equation which is similar in form to Eq. (|30D . Differences are that convergence factors have 
been inserted by hand into two of their loop integrals, a phenomenological potential has 
been used for the interaction kernel, and the perturbative expression for the gluon mass 
counterterm has been used. Their calculation also includes transverse gluon contributions 
at one loop, which is incorrect since the low momentum transverse gluon propagator is 
strongly modified by the confining interaction. As stated above, the implication of this is 
that any propagator of color nonsinglet states is IR divergent and that transverse gluons 
should only be incorporates into color singlet bound state calculations. We stress our belief 
that the interaction can, and should, be derived in a self-consistent fashion and that the use 
of perturbative renormalization is not a viable approach to nonperturbative problems. A 
related problem is seen in their Figure 2 which shows oo(q) — q approaching zero for q > 4 
GeV. This is not correct, u will in fact not approach q at large momentum, but some multiple 
of q. 

Curtis and Pennington jl4|| have noted that the gap equation derived in the Schwinger- 



Dyson formalism is not multiplicatively renormalizable. This means that it is not possible 
to exchange the A-dependence of the gap solution with a renormalization point dependence. 
As stated above, multiplicative renormalizability is not a concern in the effective field theory. 
An obvious manifestation of the difference is the appearance of noncanonical counterterms 
like the gluon mass term. In our approach the cutoff has to be kept finite, otherwise neglected 
contributions, i. e. from transverse gluons, would cause new divergences to appear. Instead 
we keep A finite and use higher dimensional local operators to account for short distance 
contributions from transverse gluons. In this way, the neglected contributions are kept finite 
(below the cutoff) and can be evaluated by expanding the Fock space used to diagonalize 
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the Hamiltonian. 

III. CONCLUSIONS 

We have presented a computational scheme for evaluating vacuum and quasiparticle 
properties of QCD. The scheme is based on the QCD Hamiltonian in Coulomb gauge, a 
field-based regularization, and the effective field theory approach to renormalization. The 
use of the Hamiltonian is advantageous because all degrees of freedom are physical and 
a transparent picture of confinement emerges - confinement is generated by the instan- 
taneous non-Abelian Coulomb interaction. We found it useful to regulate the theory by 
field-smearing because this greatly reduced the complexity of the computations. Field-based 
calculations have forced us to adopt a field-based renormalization procedure and we have 
found that the effective field theory approach is ideal for our work. 

A gap equation was derived from the QCD Hamiltonian using the Bogoliubov-Valatin 
transformation. The gap equation necessarily involves the Coulomb interaction kernel which 
should also be evaluated in the BCS vacuum. Thus self-consistent gap and effective instanta- 
neous color source interaction equations may be derived. With the aid of an approximation 
to the matrix element of the full kernel in the BCS vacuum, we have found that it is possible 
to obtain an effective potential which agrees well with the (Abelian) Coulomb and linear 
regimes of the Wilson loop measurement of the heavy source nonperturbative potential. 
Furthermore, the derived effective potential is nearly cutoff independent, indicating that a 
good low energy description has been found. 

Encouraged by this success, we look forward to repeating the calculation with the exact 
BCS interaction kernel. We fully expect this calculation to be equally successful. It will 
also be of interest to examine the consistency of our results with those derived in the quark 
sector. Since the quark and gluon vacuum sectors couple at one loop, an examination of the 
coupled gap equations should prove illuminating. Finally the efficacy of the methodology 
presented here can be tested by calculating glueball and meson spectra. 
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